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Abstract 

Some of the well-known convergence acceleration algorithms, when viewed as 
two-variable difference equations, are equivalent to discrete soliton equations. It 
is shown that the rj— algorithm is nothing but the discrete KdV equation. In ad- 
dition, one generalized version of the p— algorithm is considered to be integrable 
discretization of the cylindrical KdV equation. 
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1 Introduction 



Recently, it has been claimed that good algorithms in the field of numerical analysis 
play important roles in the theory of nonlinear integrable systems. In 1982, Symes |l] 
pointed out that one step in the QR algorithm, which is the most popular method to 
solve matrix eigenvalue problems, is equivalent to time evolution of the finite nonperiodic 
Toda equation. In 1992, Hirota, Tsujimoto, and Imai showed that the LR algorithm, 
which is another successful tool to find eigenvalues of a given matrix, is no other than 
the time-discrete Toda equation. In 1993, Papageorgiou, Grammaticos, and Ramani [|J 
showed that a well-known convergence acceleration scheme, the e— algorithm, is nothing 
but the discrete potential KdV equation. 

Our main interest in this paper is on the convergence acceleration algorithms. Let 
{S m } be a sequence of numbers which converges to Soo- In order to find by direct 
calculation, we often need a large amount of data. In such cases we transform the 
original sequence {S m } into another sequence {T m } instead of calculating directly. If 
{T m } converges to faster than {S m }, that is 

lim = 0, (1) 

we say that the transformation T : {S m } — > {T m } accelerates the convergence of the 
sequence {S m }. We now have many convergence acceleration algorithms. Among them, 
we here focus our attention on the 77—, £—, and p— algorithms [4-6]. We clarify that there 
is a strong tie between these algorithms and discrete soliton equations. 

In §2, we show that Bauer's 77— algorithm is considered to be the discrete KdV equa- 
tion in ref. M. We also look over the result by Papageorgiou et a!., the equivalence 
between Wynn's e— algorithm and the discrete potential KdV equation. In §3, we intro- 
duce a different type of algorithm, Wynn's p— algorithm. In spite of its similarity with 
the e— algorithm, it possesses noticeably different characteristics not only as a conver- 
gence accelerator but also as a discrete soliton equation. We show that the p— algorithm 
relates with the cylindrical KdV equation [|], 

u t + 6uu x + u xxx + —u = 0. (2) 

Concluding remarks are given in §4. 
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2 The rj— algorithm and the e— algorithm 



In this section we show that Bauer's //—algorithm which is one of the famous con- 
vergence acceleration algorithms, is equivalent to the discrete KdV equation. Let initial 
values 770 and rji be 



(m) 

m 
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where A is the forward difference operator given by Aa^ = a&+i — a^. Then all the other 
elements are calculated from the following recurrence relations called the rj— algorithm; 
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(rhombus rules). 
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This defines a transformation of a given series c m = r)i,m = 0, 1, 2, ... to a new series 
c' n = T]^\n = 1,2, .. . such that Y^=i c ' n converges more rapidly to the same limit S^. 
The quantities are given by the following ratios of Hankel determinants; 
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If we introduce dependent variable transformations, 
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the f]— algorithm (W) is rewritten as 



(m+1) 
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v\n -An 

which is the discrete KdV equation. Let us replace variables n and m by 

t 1 x ct 

n = -, m-- = 3 , 

e 6 lee 6 

respectively and rewrite X^ 71 ' as p + e 2 u(x — e/2, t), where e is a small parameter and c, 
p are finite constants related by 



(9) 



1 

2c=-. 



Then eq. (|||) becomes 
e 2 u(x - 1 + ce, t + e 3 ) - e 2 w(s + 1 - ce, t - e 3 



p + e 2 u(x + |,t) p + e 2 w(x - |,t) 
If we take the small limit of e, eq. ( ITU) yields the KdV equation 0, 



p 







(10) 
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at the order of e 5 . 

The discrete KdV eq. (§) is rewritten as 

{r(n + 2,m — l)r(n — 1, m) + r(n + 1, m)r(n, m — l)}r(n — 1, m + 1) 
= {t(?t, — 2,m + l)r(n + 1, m) + r(n — 1, m)r(n, m + l)}r(n + 1, m — 1), 

through a dependent variable transformation, 

x (m) _ r{n + l,m- l)r(n - 1, m) 



(13) 



(14) 



r(n, m — 1)t(ii, m) 

Subtracting r(n — 1, m+ l)r(n + 1, m — l)r(n, m) from both sides of eq. (|1"3"D, we obtain 
the bilinear form of the discrete KdV equation, 

r(n + 2, m — l)r(n — 1, m) + r(n + 1, m)r(n, m — 1) — r(n + 1, m — l)r(n, m) = 0. (15) 

It is noted that the solutions (||) and (||) are recovered by putting 
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in eq. (fH|) . 

Next, following the result by Papageorgiou et al., we briefly review the equivalence 
between the e— algorithm and the discrete potential KdV equation. The e— algorithm 
originates with Shanks and Wynn 0. Define £q and £j by 



0, e\ 
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Then all the other quantities obey the following rhombus rule called the e— algorithm; 
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et 1) )(4 m+1) -4 mj ) = i- (i9) 

According as n becomes large, e^+i converges more rapidly to Soo asm-> oo. On the 
other hand, e^n diverges as m ^ oo. This fact reminds us of the idea of the singularity 



confinement ]I0|, since a singularity at (2n, m) is confined and e^n+i converges to the 
same limit as the original sequence £i . 



It has been shown that the e— algorithm fljjp is regarded as the discrete potential 
KdV equation. The quantities s^> are also given by the following ratios of Hankel 
determinants; 
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Equation ( P0|) is called the Shanks transformation ||. Substitution of n = 1 in eq. ( pOD 
gives the well-known Aitken acceleration algorithm. 

We have so far seen that the 77— and the e— algorithms are interpreted as the discrete 
KdV and the discrete potential KdV equations, respectively. Therefore, these two algo- 
rithms are the same in their performance as convergence acceleration algorithms. This 
equivalence can also be understood from the fact [[| that the quantities and are 
related by 

(to) _ (m-1) (to) (to) _ (to) _(to-1) / 00 \ 
72n ~~ fc 2n+l fc 2n-l) '/2n+l — fc 2n+l fc 2n+l • 



3 The p— algorithm 



The p— algorithm is traced back to Thiele's rational interpolation |H |. It was first used 
as a convergence accelerator by Wynn ||. The initial values of the algorithm are given 
by 

p^ = 0, p < f ) = S m {m = 0,1,2,...), (23) 
and all the other elements fulfill the following rhombus rule; 

( P ( Z\ - pt^Wr +1] - p [ r ] ) = n. (24) 

The p— algorithm is almost the same as the e— algorithm except that "1" in the right 
hand side of eq. ( [Ill is replaced by "n" in eq. (E3|). This slight change, however, yields 
considerable differences in various aspects between these two algorithms. 

The first difference is in their performance. As one can find in ref. [P^|, the e— algorithm 



accelerates exponentially or alternatively decaying sequences, while the p— algorithm does 
rationally decaying sequences. 

The second difference is in their determinant expressions. The quantities effi are 
given by ratios of Hankel determinants, while the quantities are given by |JTT ] 



Pi m) = (-1)W%T, (25) 

n 

where [x] stands for the greatest integer less than or equal to x. Moreover, the functions 
r^" 1 ) and are expressed as the following double Casorati determinants; 
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The third difference is in the relation with discrete soliton equations. We have seen 
in the previous section that the e— algorithm is regarded as the discrete potential KdV 
equation. Before discussing the relation of the p— algorithm with soliton equations, let us 
survey the result by Papageorgiou et al. again. They considered the most general form 
of the algorithm given by 



f r M _ T ( m + 1 )w T (m+i) _ T (m)\ _ (m) 



(29) 



and applied the singularity confinement test to eq. (|29|) . As a result, when 4 m) is of the 
form, 

= f(n + m)+g(m), (30) 



eq. ( p9|) passes the test and is expected to be integrable. If we put 

f(x) = x, g(x) 



-x. 



(31) 



eq. ( p9|) gives the p— algorithm (|24|), which indicates that there is a chance for eq. (f2~4|) 
to be some discrete analogue of soliton equations. Instead of the p— algorithm Q2~3| ) itself, 
we consider the algorithm of the following form; 

(p£?i - P tX l) ){p [ n m+l) ~ P ( n ] ) = an + b(m + 1), (32) 
where a and b are constant. Employing a dependent variable transformation, 

yH=pH- p J»-i), (33) 
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from eq. (|32|) . Equation possesses a form similar to the discrete KdV eq. (||). 
However, the nonautonomous property of eq. (|34]) yields an essential difference in its 
continuous limit. Let us introduce new variables t,x defined by 
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Taking the small limit of e, we have 
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at the order of e 5 from eq. fl3~7|). Since the coefficient of u 4 is always twice as large as 
that of u/t, eq. (^) is considered as one integrable discretization of the cylindrical KdV 
equation. It is interesting to note that the p— algorithm ( p4[ ) is not exactly discretization 
of the cylindrical KdV equation. This is because we have p = in eq. ( |3T)| ) and coefficients 
of uu x and u xxx become infinite in the case of the p— algorithm. 

The third difference can be understood clearly from a viewpoint of the Hirota's for- 
malism. Employing the same dependent variable transformation as eq. (|33|), we obtain 
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from eq. (|24j) . Moreover, through the same dependent variable transformation as eq. (|T 

y(m) = F ( n + h m ~ ±)F(n - 1, m 



F(n, m — l)F(n } m) 
eq. (^) is rewritten as the following trilinear form; 

-F(n + 2,m— 1) F(n+l,m — 1) nF(n,m—l) 
F(n + l,m) F(n — l,m) 

— nF(n, m + l) F(n—l,m + l) F(n — 2,m + l) 

The functions F(n, m) and in eq. fl26[) are related by 

F(n,m) = (-ir("V( m ), 
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where a(n) satisfies 

a{n) = a(n — 2) + 



ra - 2 



(mod 2), a(0) = a(l) = 0. 
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Because of nonautonomous property of eq. ( |34"D , there is no way to derive a bilinear form 
with a single variable F(n,m) from the trilinear eq. (^T|). This fact reminds us of the 
similarity constraint of the discrete KdV equation HIBfl . It should be noted, however, that 
a pair of functions and given by eqs. (|26|) and ([27]) satisfy bilinear equations, 
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which are considered to be the Jacobi and the Pliicker identities for determinants, re- 
spectively. 



4 Concluding Remarks 

We have seen that the 77— algorithm and the e— algorithm are equivalent to the discrete 
KdV and the discrete potential KdV equations, respectively and that their performance 
as convergence acceleration algorithms is completely the same. We have also shown that 
the one generalization of the p— algorithm is considered as integrable discretization of 
the cylindrical KdV equation. The e— and the p— algorithms, despite their apparent 
similarity, possess different properties both as convergence accelerators and as discrete 
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soliton equations. The difference in performance of these two algorithms depends on 
their different determinant expressions. 

When we apply e— and p— algorithms to a convergent sequence, odd terms converge to 
the same limit as the original sequence though even terms diverge. This fact agrees with 
the idea of the singularity confinement. It is a future problem to clarify how two different 
notions, acceleration and integrability, are associated with each other. In other words, we 
should consider whether we can construct new convergence acceleration algorithms from 
the other discrete soliton equations]] and what kind of equations the other algorithms 
correspond to. The solution of these problems will shed a new light on the study of 
discrete integrable systems and numerical analysis. 
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